Precise determination of the effective-range parameters up to an arbitrary order 
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We propose a method to compute, for a given potential model, an arbitrary coefficient of the 
effective-range function expanded as a power series in energy. The method is based on a set of 
recurrence relations at low energy, that allows a compact and general description to any order 
in energy for neutral and charged cases. By using the Lagrange mesh technique to compute the 
R-matrix at zero energy, this proposal permits to compute, with a very good precision, the effective- 
range parameters. We use a potential model for some nuclear systems to illustrate the effectiveness 
of this method and to discuss its numerical limitations. 



I. INTRODUCTION 

Since the early days of nuclear physics, the study of 
scattering processes at low energies has been based on 
the partial-wave analysis of cross sections in terms of 
phase-shifts <5/, where / is the orbital quantum number. 
Among the first systems studied, the description of the 
s-wave elastic scattering for the neutron-proton system 
took a high interest [l|, |2j. For this system in particular, 
a good approximation of the function fccot((5o) (where 
k is the wave number in the center-of-mass frame) as 
a linear function of k 2 led to study the Effective-Range 
Function (ERF) and its expansion in Maclaurin series, 
better known as the Effective- Range Expansion (ERE). 
Since then, an ERF has been established for arbitrary 
partial waves, both in the neutral and charged cases [3|. 
The analyticity properties of this function have also been 
demonstrated, making of EREs an important tool for the 
parametrization of low-energy scattering data. 

The coefficients of the ERE, which are directly related 
to the so-called effective-range parameters, have become 
a field of study in themselves |j-[6( . These parameters, in- 
cluding the well-known scattering length, effective range 
and shape parameter, can be interpreted in terms of 
physical properties of the system, in particular regard- 
ing scattering states at low energies and weakly bound 
states Q- They can be considered as parameters fitting 
these properties, hence providing a useful analysis tool 
for experimental data. For instance, Ref. [2| shows a way 
to obtain the first parameters of the ERE for the neutron- 
proton and proton-proton systems from experiment. Of 
course, in many cases, the possibility of determining the 
effective-range parameters via the experimental data is 
open but the precision will depend on how close we can 
go to zero energy and how many parameters we wish to 
determine. 

The ERE coefficients actually carry information about 
the interaction between the colliding particles; they can 
thus be considered as an intermediate step between ex- 
periment and models. In an inverse-problem perspective, 
they can be extracted from scattering data and then used 
as an input for an inversion method aimed at deriving an 
interaction model. Conversely, ERE coefficients are cal- 
culated for a given interaction by solving the Schrodingcr 



equation at low-energy scattering, hence providing an ef- 
ficient way of comparing the predictions from this model 
with experimental data. A efficient method to obtain 
the first three effective-range parameters is proposed in 
Rcfs. [a Q. In particular, the latter shows that using 
the R-matrix [Hj, [llj and the Lagrange mesh technique 
[la [l3| , the prediction of the scattering length, the ef- 
fective range and the shape parameters is accurate for 
neutral and charged cases. The aim of the present paper 
is to extend that proposal to calculate any parameter of 
the ERE in both cases, for a given interaction model. 

In Sec. |lll we introduce the ERF of a two-body sys- 
tem via the phase shift both for the neutral and charged 
case. Then, we define the effective-range parameters and 
present the method to compute any of them. In Sec. IIIII 
we briefly describe the scheme to compute derivatives at 
zero energy of the R-matrix based on the Lagrange mesh 
technique. In order to test the method, three neutral and 
three charged cases are analyzed in Sec. IIVI Conclusions 
and perspectives are presented in Sec. [Vj 



II. EFFECTIVE-RANGE FUNCTION AT LOW 
ENERGIES 

Let us consider the scattering of two particles with 
relative coordinate r and reduced mass /i. Let us also 
suppose that part of the interaction between the two 
particles decreases fast enough in such a way that the 
effective potential can be described in two regions as the 
contribution of: a short-range potential (Vn), a centrifu- 
gal potential and a Coulomb potential (Vc) for r < a 
and a centrifugal and a Coulomb potentials for r > a. 
Here a is considered as the range of the short-range po- 
tential but in practice it could be taken larger or equal 
to the minimal distance between the two particles when 
the short-range potential is neglected in comparison with 
the rest of the interaction. In this scenario, for the Zth 
partial wave, and if ^i(k, r) = - fc ^. is the radial wave 
function, the radial Schrodinger equation can be written 
as 



tf_£_ 

2n dr 2 



h 2 l(i + i) 



+ V(r) - E 



u,(fc,r)=0, (1) 



where ui(k, r) is known as the modified radial wave func- 
tion, E = |-fc 2 is the energy in the center of mass frame 



and 



V(r) 



Vv(r) + V c (r) 
Vc(r), 



if r < a, 
if r > a. 



(2) 



Let us introduce the following notations and conven- 
tions: 

• Hi; = Hn(k,a) is the Ricatti-Bessel (regular 
Coulomb) function at r — a for the neutral (charged) 
case. 

• H 2 i = H2i{k,a) is the Ricatti- Neumann (irregular 
Coulomb) function at r = a for the neutral (charged) 
case. 

• Prime represents the partial derivative with respect 
to r at r = a (e.g., F^ = dH u (k,r)/dr\ r=a ). 



R 



matrix [1 



dln( Ul {k,r)) 
dr 

a. 



is the single channel R- 



• For the charged case, if the interacting particles have 



electric charge Z\e and ^e, then ajy 

h 2 1 



4-n-ep 



fj ZiZ 2 o 2 



the nuclear Bohr radius, En = 5-t- is the nuclear 



Rydberg energy and 



1 



a^k 



(3) 



is the Sommerfeld parameter. 



With the previous notations and conventions, the phase- 
shift Si J4| is obtained through 



tan(<5; 



ill j — aRiH' u 
aRiH 2l 



H- 



2/ 



(4) 



Equation Q gives a good tool to compute the phase- 
shifts and later the ERF in a wide range of energy. Un- 
fortunately Eq. ((31) is not useful at low energies, in partic- 
ular at the limit when E tends towards zero, because the 
functions Hu and H21 and their radial derivatives van- 
ish, giving an indeterminate value for the phase-shifts 
and then for the ERF. However, this limit is very impor- 
tant to determine the coefficients of the ERE, since by 
definition they are calculated at zero energy. 

In order to overcome the indeterminacy of Si at E = 0, 
a convenient re- normalization of the Hu and H21 func- 
tions is proposed in Ref. [9j. These renormalizcd func- 
tions remain finite at zero energy and are expressed as 



W-n — 



H-21 



k Hu, (neutral case), 

k~ x / 2 exp(Trrj)Hi h (charged case), 



k l H 2h 
ffc- 1 /2 C xp(-7rr / ) J ff 2 ,, 



(neutral case), 
(charged case). 



(5) 
(6) 



Using definitions ([5]) and (J6)) in Eq. (|4]) one can easily 
obtain the expressions for the phase-shits 



tan(<y 



k 2l+1 D u 



(neutral case), 



§ exp(— 2m])Di, (charged case), 



with 



Di = - 



Hu - aRiH' u 
H21 - aRiW 



(7) 



(8) 



2/ 



For simplicity, we have omitted the dependence on energy 
inEqs. © and ©. 

By using the previous notation the ERF, Ki, can be 
written as [3[ 



K, 



d; 



2D; 



h] 



V 2 a 2l + 1 

where the functions wi and h read 



(neutral case), 
(charged case) 



wi 



1 

n 

3=0 



i/}(ir]) — ln(iry) 



1 + ^ 
9 

1 



2i v e 2*v - 1 ' 



(9) 

(10) 
(11) 



with tp the digamma function [14J . 

It can be proven that Ki is analytical at zero energy, 
which allows to write this function as a power series in 
energy, i.e., 



Ki = J2 c ^ Er 



(12) 



71=0 



Equation (fT2")) is useful to define the effective-range pa- 
rameters. Here we adopt the standard convention for 
the first three efective-range parameters: the scattering 
length (ai), the effective range (rj) and the shape param- 
eter (Pi), i.e., expanding Eq. (fT2")) one finds 



Ki 



1 

ai 

■£■ 

n— 3 



__ 
2 

cl,n 



2/' 

w 

2//. 



E 



Pirf 



E 2 



E' 



(13) 



/ 2 \ n 
In Eq. ([T3| we have defined Qi :U — ( |- ) Q,« as the 

effective-range parameter of order n, for n > 3. 

Our goal is to compute the coefficient c; >n in Eq. (TT2"]) 
and then to calculate the effective-range parameters in 
Eq. (fl"3"l) . The coefficient c; >n is easily found by computing 
the nth derivative of Ki at zero energy 



Q. 



f(^)o' 



(») 



(14) 



Hereafter, we use the superscript (n) to designate the nth 
partial derivative with respect to energy and the subscript 



to designate functions at E = 0. The reader should not 
confuse the prime with the superscript (1). 

As we have shown in Eq. (J9J , the explicit form of Ki 
depends on whether the Coulomb interaction is present 
in the system or not. For this reason we will study both 
cases separately in the next two subsections. 



A. ci t n f° r the neutral case 

In the neutral case c^ n depends on _D ; _ only, then 
using Eq. (|8]) one can rewrite Eq. (fl4|) as 



Q. 



L / H21 — aRiH 2 i 



n\ \ Hn — aR{M/ Xi 



(15) 



// 



i\(") 



(°rX 



-U(n) 



The term (Df ) can be expanded in terms of the 
derivatives of Hit, W u and Ri at E = 0. To do this, 
let us first define D^ 1 as the following 



(H 2l t ) = (^y(2l-2n-l)\\a- l+2 * 



(22) 



Thus, there only remains to find a general expression 
for (Ri)q to complete all requirements to compute the 
coefficient ci :U . An efficient method to get the derivatives 
of the R-matrix at zero energy is described in Ref. Q ■ We 
will briefly explain it in Sec. 



B. c; n for the charged case 

By using Eqs. ((9]) and ([T3]l . the coefficient ci >n for the 
charged case reads 



Q. 



n\ Z! 2 < +1 



>-i\(«) 



2(vnD r X +(ww 



(n) 



(23) 



Note that one can use the same procedure as the one used 
for the neutral case to split the derivatives of the products 

wiD^ and w;/i. With this idea the term (wiD^ ) can 
be written as 



D,- 1 = -^ 



A2 
'Ai, : 



(16) 



with 



A it = H %1 - aRM'n, for * G {1, 2}. (17) 

Using the product rule for derivatives, one gets 

-l\( n ) _ V^ I n \ ( A \(n-m) /* -l\(m) 



D 



m= 



with ( ra ) the binomial coefficient. In the same way, the 
terms (A*j)q and (A i; _1 ) can be written as 



(AuP = (n tl )i n) - «E ( n ;wt~ i HK)P (19) 

j=0 ^ J/ 

and 



n n—ra , -i \ / -\ 

-i\(n) _ sr^ v^ l n ~ l \ (n-m\ ^ A .( m ) 



(**tw =-y,j: 



ni — 1 j—0 



111 — 1 



(A 



iUa 



-l\(n-m-j) 



(A," 1 ) 



,M 



(Ar 1 )^, (20) 



for n > 0. 

The four previous equations, together with Eq. (1151) . 
allow us to compute recursively any coefficient ci jTl if 

one knows (Hu)q , (^i;)o an d (Ri)o ■ The value of 
(H'n) is easily obtained from (Hu) since the deriva- 
tives over r and 2£ can be interchanged. The functions 
(H u i n) and (« 2I )W rea d @ 



(^iz)o' 



(n) 



fi\« a 



i+2n + l 



ft 2 / (2; + 2n+l)!!' 



(21) 



(wj-D 



-i\( n ) 
1 )o 



E 



M 



(»-m) (A -1 } (-) (24) 



and using the power series expansion of the function h 

h 



^ 1 E IE 2 IE 3 IE 4 

> h n E = — — 



A' 



120^2 252 E\ T 240^1 



J JV 



JV 



1 E 5 691 £ 6 1 E r 
132 £^ + 32760 E% + 12 ^ 
3617 £ 8 43867 £ 9 
8160 £^~ + 14364 £^ + '" ' 



J N 



(25) 



the term (wih) n reads 



(n) 



(«*4 n) 



^-"' (n — m)! 

m=0 v ' 



MS 



(n— m) j 



(26) 



The term (Z?j x ) ™ in Eq. (|2"4"jl can be computed by 
using Eqs. ([18 ]) - ([2"0"]l . For n > the term (w;)o reads 

n-l 



M{, 
with 



(«) 



I ill / \(ra — 1— m) 

V^/ lV n (n-l)! (w;)^> 

Z^f." 1 ^ 7 ; - ^m iN %™< 



(n — 1 — m) ! 



3=1 



pm+l 



(27) 
(28) 



v(0) 



The value (wi) = 1 is easily calculated from Eq. ([10]) . 
In Eq. (|25[) we have written the first ten terms explic- 
itly, which will be used in Sec. IIVI to compute the first 
ten effective-range parameters. Note that ho = 0. 



At this point we have all the requirements to compute 
Ci tn except for general expressions of (Ri)q , (H'h)q and 
{Hh)q which arc required to calculate Eq. (Tl9l) . Simi- 
larly to the neutral case, the term (Ri)q will be obtained 
by using the Lagrange mesh technique described in Sec. 
IIII1 and the term {H' u )q will be easily computed if one 
knows CHu)o ■ Therefore, we can focus on finding a 

general expression for (Hu)q ■ 

Let us start with the regular Coulomb function -Hi;. 
For 77 > this function can be expanded as 



H u 



ii'i 



2-q e 2m > - 1 



2/ 



E * 

m=2l+l 



[2) T m(x), 



(29) 

where we have started from its expansion in terms of 
Besscl-Clifford functions [Tj]. The argument of the first 
modified Bcsscl function I m is the dimensionlcss variable 

— . The term b m satisfies the recurrence relation 



4r) 2 (m - 2l)b m+ i + mb m -i + b m -i 



0. 



(30) 



for m > 2/ + 1 and with bii = b^i+2 = 0, ^22+1 = 1- 

By using Eqs. (f2U| in Eq. (|5|) one finds the compact 
expression for the renormalized Coulomb function 



Hu = Vnr ^2 d m (f> m (x,l), 

m=2l+l 



(31) 



where the terms d n 



d m (r],l) and <j) m (x,l) read 



is given by the term b my /wi. Therefore, from Eq. ([31 
one obtains 



CHu) 







E (bmVm)o' )( t> m (x,l). (34) 



ro=2Z+l 



Using the product rule for derivatives and taking into 
account that b m is a polynomial in energy, Eq. (|34[) can 
be written as 



n 3(ra-j) 

(Hup = v^y: e 

3=0 m=2(n-j) 



y (^)« 



x(b 



m+2l + l 



(n-j) ( x 
2 



'0 



=r) I, 



m+2l+l 



(■r) 



(35) 



Note that the sum over m runs from 2(n — j) up to 3(n — 
j). This is because (& m +2M-i)o = if m < 2m or 
m > 3m, as one can prove from Eq. (|30l) . 

The terms (v^) and (6 m+ i)o can be calculated 
by iterating 



(V550 



(„) _ (n- 1)! 



(i 



2B5 



vw 7=0 
with iVj j defined in Eq 



Ec- 1 



,(n-l-i) 



(n-l-i)! ( ' J 



^ij. ( 36 ) 



and 



n \(n— 1) 1 /? \(n— 1) 

(n) _ M "l(0m-l)o +(Om-2)o 



(&m+l)o = -n 



4(m - 20 -E. 



N 



(37) 



«'; 



1 - er 2 ^ 



(32) 



for 771 > 1 and 77 > 0. 

For i € {1,2}, the simple and useful relation 



and 



m (x,O= (2) 7 ™(z)- ( 33 ) 

Because of the fast decrease of the term e -27 ^ when the 
energy goes to zero, one can check from Eq. (|32[) that the 
only contribution to the derivatives of d m at zero energy 

I 



/ s(n) 



(Hii)a 



1 
27 



w x ,W|? 



(Hii) JrX 



<:).r 



(38) 



helps us to compute (H'h)q easily. 

If i = 1 and using Eq. (|35|) one can compute the term 

a( ^""' in Eq. J35J) as 



d(Hut ] 

dx 



n 3(n-j) 

v^E E 

j=0 m=2(n-j) 



(Vwi)o (bm+2l+l)c 



m , . d/ m+ 2i+i(a;) 
-I m+2 * +1 (z) + 



(39) 



We have discussed the steps to obtain (Hh)q and 

(H'il)o if i — 1- We can now move on to the case i = 2. 
For this case, it is necessary to clarify that the expan- 
sion in Bessel-Clifford functions of the irregular Coulomb 
wave function H 21 [14| is possible only when 77 — > 00. This 



I 

limit is equivalent to E — > (see Eq. ((3])) which is our 
case of interest. Therefore, mimicking the procedure fol- 
lowed previously for (Hu)q and {H' u )q , and defining 
I n as the second modified Bcsscl function, we find 



m^p 



n 3(n-j) 



dx 



(») 



j=0 m=2(n-j) 
n 3(n—j) 



j=0 m=2(n—j) 



■[".) (V^i)^ (b m+ 2i + it- j) (1)™ Wm+i( 



Uj) 



(-1)™( ? )(V^)^(6 Tn+2;+1 )^ 



(n-j) f •£ 

2 



— lm+2l+l(X) 
X 



hl(z) 


) 


(40) 


dI m +2Z+l( a; ) N \ 


■ (41) 




dx / 



III. LAGRANGE MESH TECHNIQUE TO 

COMPUTE R-MATRIX DERIVATIVES AT ZERO 

ENERGY 

In order to understand some properties of two-body 
nuclear systems, it is convenient to take advantage of the 
fact that the nuclear interaction is a short-range inter- 
action. This fact allows us to divide the configuration 
space in two regions as it is shown in Eqs. (fl]) and ©. 
In the external region (r > a) the Schrodinger equation 
is solved by analytic procedures. The challenge is to find 
the wave function in the internal region (r < a) together 
with the correct matching between the external and in- 
ternal wave functions at r = a. It is equivalent to find 
the R-matrix for a single-channel case. In general, using 
analytic methods it is not possible to achieve this chal- 
lenge (at least nowadays) and for this reason one needs to 
find a good approximation of the wave function, which is 
usually done numerically. In the case of the R-matrix Ri , 
Ref. ll| shows how to achieve this challenge by approxi- 
mating the wave function as a superposition of N linearly 
independent functions, /«(?"), of a basis. The result for a 
two-body system with reduced mass (i and energy E is 



A' 



Ri 



2fia 2 



J2 fn(a)[(C - Eiy l ] nm f m (a), (42) 



where the elements of the N x N matrix / and C are 
respectively the Kroncckcr delta S nm and 



' a 



U(r)[Ti + V{r) + £]f m (r)dr, 



with 



Ti = 

2(i 



dr 2 



1(1 + 1) 



(43) 



(44) 



the kinetic energy operator of the partial wave I, and 

(45) 



2^ (r_o) *' 



the Bloch operator [15J . The latter appears for two main 
reasons. First, "to keep" the hamiltonian hermitian in 
the internal region and second, to match the internal and 
external wave functions correctly. 

The good estimation of G nm in Eq. (|4"3")l will depend 
on the basis chosen and how the integral is made. As 



the internal space is r £ [ 0, a] , a good basis is the shifted 
Lagrange-Legendre basis (ll[. This basis allows to com- 
pute the integral in Eq. f|43[) easily if it is approximated 
by the shifted Gauss-Lcgendre quadrature. With these 
assumptions / n (a) and C nm read [9l. fill] 



fn(a) 



(-1)" 



for n£ {1,2,--- ,N}, (46) 



L ii 



\Jax n (\ - x n )' 

( h 2 (4A 2 +4A+3)a:„(l-a:„)-6a:„ + l 
2ft 3a 2 xl(l-x„) 2 

+ i)l l '^xY + v ( ax n), for n = m, 

h 2 (_!)"+"■ 



2,1 a 2 y/ ' x n x m (l-x n ){l-x m ) 

r2 i AT i 1 i x n -\-x rn —2x n x rl 



(47) 



N 2 + N+1 + 



(x n —Xm)' 2 

, for n =/= m, 



with x n the abscissas of a shifted Gauss-Lcgendre quadra- 
ture in the [0,1] interval. 

Finally, the previous approximations allow to calculate 
the jth derivative of the R-matrix at zero energy as 



(Rii j) 



2/i a 



A 

E 

n,m— 1 



fn(a){(C-') 



-iy+ii 



J m (a). (48) 



IV. APPLICATIONS 

In this section we discuss the effectiveness of the 
method derived in Sec. |TT] to calculate the first ten 
cffcctivc-rangc parameters via the coefficients a. n in Eqs. 
(fT5]) and (|2"3"|) . To do that, we use the technique shown 
in Sec. IIIII to compute the R-matrix derivatives at zero 
energy. We analyze the following nuclear systems: 

Neutral systems: 

• SW (Square Well): s,p, d-waves, |- = 1 MeV fm 2 and 
square well potential of width 3 fm and depth 10 MeV. 

• W-S (Woods-Saxon): s,p,d-w&ves, f£ = 1 MeV fm 2 
and Woods-Saxon potential 1+ex p iv^m/ d i > "vvith R = 3 
fm and V = 10[1 +cxp(-R/d)] MeV. The diffuscness 
d will be taken as a variable parameter. 

• n-p (neutron-proton): s- wave and Bargmann potential 



i 

fm 



cxp(-2/3r) 
[l+6cxp(-2,3r)] 2 2/i> 

1 and 8 = 0.81 fm" 1 Ki 



with b 



P-a 

/3 + Q' 



0.04 



Charged systems: 



a+ He: s-wavc, point-sphere 
with R c = 3.248 fm and 

r.2 



7(2.477 im) 2 



Coulomb 
gaussian 
MeV 1. 



— 67.67 exp 

16 0+p: p-wave, point-sphere 
with Rc = 3.553 fm and 
-36 exp [-r 2 /(3.553 fm) 2 ] MeV [f. 

12 C+a: d-wdve, erf-Coulomb 



Coulomb 
gaussian 



potential 
potential 

potential 
potential 

potential 



Z X Z^ 



crf( 



47reor ^ iA V 2.5 fm 

and gaussian potential 
MeV 0. 



where erf is the error function [14j 
112.3319 exp [-r 2 /(2.8 fm) 2 ] 



All the previous systems are chosen in order to explore 
the numerical merits and numerical limitations of our 
method as much as possible. We choose the SW and n- 
p systems because they allow us to compare our results 
with exact results. The W-S system is chosen to evaluate 
the behavior of the effective-range parameters when a 
weakly bound state, virtual state or resonance at low 
energy is present. The last three systems were chosen 
to evaluate the charged case and because, together with 
n-p, all of them are closer to real nuclear cases. 

Before showing results, let us specify the calculation 
conditions. All masses are atomic masses and are taken 
from Ref. [17| except for systems SW and W-S. Units of 
the effective-range parameters and theirs absolute errors, 
e, will be given in fm" where v is shown in Table |U 



Table I. Units for the effective- 
error in fm". 


-range 


parameters and absolute 




ai, e(ai) 


ri, e(n) 


Pi. 


<Pi) 


Ql.n, s{Ql, 


n) 


V 


21 + 1 


-21 + 1 




41 


-21 + 2n- 


1 



In order to estimate the effective-range parameters us- 
ing the method described in Sec. HH it is necessary to 
compute expressions that involve many operations that 
have to be carried out numerically. This fact leads to 
some numerical limitations in the accuracy of the esti- 
mations that will depend on the characteristics of the 
software and hardware used. Here we develop such a nu- 
merical evaluation using double precision in FORTRAN 
language on a regular desktop computer. With these con- 
ditions one expects that the values estimated match the 
exact result up to 15 digits in the best case. 



A. Neutral cases 

Exact values of the first ten effective range parameters 
for the SW and n-p systems are displayed in Table QI] 

Let us discuss some relevant features of the Qi t „ pa- 
rameters displayed in Table [H] from the mathematical 
and numerical point of view. First, for the s, p and d 
waves of the SW system, the absolute value of each pa- 
rameter decreases with its order for all cases except for 



Table II. Exact values of the first ten effective-range parame- 
ters for the SW and n-p systems. All values are written with 
twelve digits except Po and Qo,n for the n-p system. Units 
can be consulted in Table [I] 

SW 



s-wave 



p-wave 



d-wave 



<7 

Pi 

Ql,3 

QlA 

Qifi 

Ql,6 

Ql,7 
Ql,8 
Ql.9 



2.98035125196 
1.95321791835 
-0.07224549975 
0.41010095011 
0.32190586543 
0.25278100437 
0.19808608111 
0.15503303001 
0.12127431544 
0.09484976310 



54.5044450592 
-1.01095733502 
0.69171350459 
0.21472158944 
0.10883657253 
0.04415761104 
0.02925502308 
0.00791996459 
0.00948504977 
-0.00012111901 



4.32199541189 
-1.04292323145 
-0.87150574727 
0.35035876175 
0.01058724882 
-0.01015831203 
-0.00625258947 
-0.00253456328 
-0.00082189556 
-0.00021419743 



n-p 



a 

''o 
Po 

Qo.n 



-23.7654320988 

2.59740259740 



0, for n > 3 



Qi,8- This behavior guarantees a good convergence at 
low energies for the ERE, and then, a good description 
of the ERF and phase-shifts with a small set of effective- 
range parameters. For the n-p system this effect is rather 
evident because the ERE converges for all energies. Sec- 
ond, for the s-wavc of the SW system, all Qi. n parame- 
ters in Table [TT| have the same sign and almost the same 
order of magnitude. This fact restricts the range of con- 
vergence of the ERE but the numerical estimation for 
the parameters should be equally good or bad for all, if 
the numerical precision can handle correctly all calcula- 
tions of the method, i.e., the results should be reliable 
except for cases when subtractions or additions of big 
numbers (from the numerical point of view) are involved 
in the calculation. In contrast, this explanation could 
play a relevant role in the calculation of Q;„ for the n-p 
system, where the numerical calculation should compute 
zero from subtractions and additions that maybe involve 
big numbers. Third, for the d-wave of SW, the parame- 
ters do not have the same sign and the order of magnitude 
decreases every second order. This increases the range of 
convergence of the ERE for the d-wave. 

Let us now show and analyze some results for SW by 
using the method in Sec. HU For all cases (s,p, <i-waves) 
there is only one acceptable value for the channel radius, 
a = 3 fm. For a < 3 fm, part of the potential is neglected 
and therefore the R-matrix does not give physical results. 
For a > 3 fm, the Lagrange mesh technique, and particu- 
larly the Gauss quadrature, cannot handle correctly the 
discontinuity in the potential at r = 3 fm and therefore 
the R-matrix is wrongly estimated. Thus, choosing a = 3 
fm, Table ITO1 displays the absolute error for the first ten 
effective-range parameters for s,p, d- waves of SW. 

Tabic Mil highlights a good estimation of the effective- 
range parameters for N = 10 and very good for N > 20. 



Table III. Absolute errors of the first ten effective-range pa- 
rameters for s,p, d- waves of SW. The absolute error is written 



in normalized scientific notation as e 



x 10 , where 



is the significand and b is the order of magnitude, 
report b only. Units are shown in Table [J 
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This is expected because the square well potential has 
a polynomial form in the region r < 3 fm, and then, 
the Lagrange mesh technique and Gauss quadrature give 
better approximations for the R-matrix than for those 
potentials with a different dependency on r. 

A relevant feature in Table IIIII is the remarkable in- 
creasing of e(Qi.s) in comparison with those of lower or- 
der. This effect could be expected because Qi.g is almost 
two orders of magnitude smaller than the parameters of 
lower order, which implies that its numerical calculation 
requires higher computational precision. On the other 
hand, the ninth order correction in the ERE does not 
play a relevant role, which can be corroborated by com- 
paring the exact phase-shift and its calculation using the 
ERE up to eighth and ninth order. This is another way to 
check whether an effective-range parameter could require 
higher numerical precision than others of lower order. 

Note that in Table [TT] the absolute values of the 
effective-range parameters for the p-wave of SW look sim- 
ilar to those for the <i-wave but the estimation is worse 
for the p-wave as Table IIIII shows. It leads us to think 
that there is an extra feature that makes poorer estima- 
tion for the p-wave case (especially for high order pa- 
rameters). This feature is related with the presence of 
a weakly bound state, virtual state or resonance at low 
energy. We shall discuss their effects later in this section. 

At this point we have shown that our method in combi- 
nation with the Lagrange mesh technique works very well 
for simple potentials. Now we wish to move on to more 
sophisticated and realistic shapes for the short-range po- 
tential. Let us first consider the W-S system. In this case 
we do not know the exact values of the effective-range pa- 
rameters to compare with our estimations, therefore in 
order to have a check point, we shall start with potentials 
close to the SW potential as Fig. Q] shows. 

The W-S potential does not vanish for a finite distance, 
which implies that there is not a precise value for a. How- 
ever, the Gauss quadrature should work correctly from 
the mathematical point of view because this potential 
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Figure 1. Short-range potential Vjv plus centrifugal potential 
Vi for the SW (d — 0) and W-S systems. Three examples are 
shown for the W-S system. Labels give the diffuseness in fm. 



has no discontinuity for d > 0. 

As we discussed at the beginning, the restriction on a 
is that the short-range potential could be neglected for 
r > a. Thus, in order to obtain acceptable results, we 
scan the channel radius from a minimal value a obtained 
from | Vat (a) | = | Vjv(0)/100|. We also scan the diffuseness 
in the interval d < 0.6 fm and choose N — 30 for all cases. 

In general, we observe a smooth and slow change in 
the effective-range parameters when the diffuseness in- 
creases. Even more, they reach a stable value for chan- 
nel radii close to a. An exception to this behavior occurs 
when a weakly bound state, virtual state or resonance at 
low energy is present. It is shown in Fig. [2] where we have 
chosen the scattering length as an example. 
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Figure 2. Prediction of ai for the W-S system by choosing 
different channel radii. For all cases TV = 30 except for the 
diamonds in the middle-left plot where N = 40. Labels give 
the diffuseness in fm. The lines are included to guide the eye. 



Remember that the choice of the channel radius implies 
a truncation of the potential range for potentials that do 
not vanish at a finite distance. This fact in practice, leads 
to a kind of domino effect, i.e., channel radius — > narrower 
potential — > slight rise of the energy levels — > chance to 
modify a weakly bound state making it a virtual state or 
resonance at low energy — ¥ incorrect estimation of the 
effective-range parameters. Thus, in order to prevent 
wrong predictions, the channel radius should be large 
enough to keep the spectrum unchanged, especially when 
a weakly bound state is present. 

In order to discuss the impact of virtual states and 
resonances at low energy and to contrast their effects 
with those given by weakly bound states, we show the 
energy levels of the W-S system for different values of d 
in Fig. El 
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Figure 3. Scheme of energy levels for the systems W-S and 
SW (d — 0). The diffuseness is shown in parenthesis and is 
given in fm. Solid line: I — 0. Dashed line: 1 = 1. Dashed- 
dotted line: I = 2. 

From Fig. [3J the most evident cases when a weakly 
bound state is present are d = 0.1 fm for the p-wave and 
d = 0.6 fm for the s-wave (we skip the case d = for the 
s-wave, which corresponds to the SW system). Looking 
at Fig.[5]for these two cases, one infers that weakly bound 
states do not introduce new restrictions for the channel 
radius to those based on the domino effect discussed pre- 
viously, but may demand a better approximation of the 
wave function, as it is illustrated by both curves in the 
case d = 0.1 and p-wave in Fig. [2] 

On the other hand, comparing the s-wave cross sec- 
tions at zero energy of the W-S system for d = 0.4 fm 
(cto(0) « 0.1 b), d = 0.5 fm (a (0) ra 10 4 b) and d = 0.6 
fm (i7o(0) « 10 b) one sees a remarkable increase of this 
for the case d = 0.5 fm. Taking into account that there 
is no weakly bound state for d = 0.5 fm as Fig. [3] shows, 
then this increase is due to the presence of a virtual state 
at low energy. The effect of that virtual state on the es- 
timate of ao is shown on the top-right panel of Fig. [2] 
where the channel radius must be larger than the range 
that could be chosen by simple inspection of Fig. [TJ The 
same effect is seen for the <i-wave and a diffuseness of 0.6 



fm. In this case, the d-wave cross section shows a reso- 
nance around 66 keV which leads to large values of a as 
the bottom-right panel of Fig. [5] shows. 

Note that for the cases where there is neither a weakly 
bound state, nor a virtual state or resonance at low en- 
ergy, Fig. [5] shows that the estimation of ai (and in gen- 
eral the effective-range parameters) is more stable. This 
leads us to summarize the constraints on the channel ra- 
dius to correctly compute effective-range parameters as: 
1) the presence of a weakly bound state can demand a 
channel radius larger than the typical range of the po- 
tential (see for instance the definition of a) and 2) if a 
virtual state or resonance at low energy is present the 
channel radius must be much larger. 

We wish to clarify that the previous constraints are 
based on short-range potentials which arc not null but 
negligible for r > a. For potentials strictly vanishing be- 
yond a finite distance, a very good estimation is expected 
when a is equal to the potential range (see for instance 
Table liLTl and Fig. [3] for the p-wave of the SW system). 

For the n-p system we also know the exact value of 
each effective-range parameter. It makes this system a 
good candidate to test our method. In Fig. @]we show 
the potential shape for the n-p system and for the three 
charged cases which will be analyzed later. 


-10 



- 



d 



-20 - 
-30 - 

-40 

-50 

-60 

-70 




■ | ' 'f O+p'^^-- 1 




1 


1 


..y-' 


'"* ,.-" : '"'" 


- 


i A p / 






- 


J Q+ 3 He/ 






_ 


V / 2 C+a 






- 


— '' | "" | | 




, 





12 3 4 

r (fm) 



Figure 4. Short-range potential Viv plus centrifugal potential 
Vi for systems n-p (s-wave), a+ He (s-wave), O+p (p-wave) 
and 12 C+a (d-wave). 

Tabic IIVI displays the absolute error of the first ten 
effective-range parameters for the n-p system. Note that, 
for N G {20,30} and a G {10, 12, 15} fm, the scattering 
length, the effective range and the shape parameter are 
in perfect agreement or slightly better estimated than 
the results of Ref. |{|. The estimation for the first six 
effective-range parameters is very good if a = 17 fm and 
N = 40. However, the precision for the last four parame- 
ters is very bad in general, even increasing the Lagrange 
mesh points. The possibility of increasing a does not give 
us better results as one sees in Table HVl for a = 19 fm. 

The reader could think that a = 17 fm is very large 
following Fig. SI which shows that 4 to 6 fm seems good 
enough. This increase in the channel radius is expected 



Table IV. Absolute errors of the first ten effective-range parameters for the n-p system. The absolute error is written in 
normalized scientific notation as s = e s i g x 10 6 , where e B - lg is the significand and b is the order of magnitude. Here we report b 
only. Units are shown in Table U 
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because the Bargmann potential used for the n-p system 
provides a virtual state at —66 keV. This makes the n-p 
system a good example where the choice of the channel 
radius should take into account the presence of a virtual 
state at low energy. 

The wrong calculation for the last four parameters 
shown in Table lFvH are explained from the numerical point 
of view. In those cases, double precision is not enough 
to keep numbers with a large set of digits, which are 
needed to compute small quantities from subtractions of 
big ones. It can be illustrated by comparing interme- 
diate results from the numerical calculation with those 
obtained analytically as Fig. [5] shows for the R-matrix 
and its derivatives. 



a: 




Figure 5. R-matrix derivatives at zero energy obtained from 
the scattering s-wave function for the n-p system (solid line) 
and from numerical calculation of Eq. (|48[) (square points). 
Black (white) points correspond to N = 40 (N = 20). Four 
different derivative-orders are shown as examples. 



In this figure (Ri)q is computed analytically by us- 
ing the scattering wave function for a Bargmann poten- 
tial |5|, [l6j and numerically by using the Lagrange-mesh 
technique. 

Figure [5] highlights that the R-matrix at zero energy 
(j = 0) is very well computed numerically for different 
Lagrange meshes. It agrees with the results for ao shown 
in Table IIV1 where one can deduce that the difference 
between the value of ao for TV = 20 and N = 40 is given 
beyond the fourth digit (remember, the exact value is 
ciq = —23.765 • • • fm). Similarly for the R-matrix deriva- 
tives at zero energy, the bottom three panels of Fig. [5] 
show that the numerical results are in a good agreement 
with the exact values. 

As the R-matrix derivatives are well computed and 
the numerical calculation of the (Hh)q is accurate, 
the wrong estimation of the effective-range parameters 
should come from computational processes. In order 
to explain it, we wish to emphasize that these param- 
eters are obtained in three steps: First, compute ("Hiz) 
and (Ri)i j \ Second, calculate Eqs. Ifity , (|2TJj) and (fig]) . 
Third, determine the effective-range parameters via Eq. 
(fl~5|) . Let us detail the third step by expandig Eq. (fT5")) 
as 



Cl,n 



»!(Aii)§ 



(7^1 z)o™ C^2/)o - (%2;)o™ 0#l/)o 



+ a(R l )$ l \H' 2 M'Hii)o 



a(Ri)^(H' u ) (^2i)o- 



(49) 



Equation (|49[) shows that numerical limitations can oc- 
cur if large values of t\ , t%, t^ or T4 (or any of the remain- 
ing terms) appear for a given numerical precision. For 
instance, by using double precision and choosing a = 17 
fm for the n-p system, one finds that to calculate Co,g, 
which exact value is zero, the terms ti«2x 10~ 3 M f "y U , 
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To 



4xl0- 2 ^aiidT3 



M f ^a do not contribute to 

5X1 TTTTrEr 



the numerical imprecision as the term T4 ~ 
does (the reader can check these values by using Eqs. 
(f2"Tj) . (J22) and Fig. [5]). Comparing t 3 and t 4 one sees 
that double precision is not enough to compute correctly 
Co. 9 and then Qo,9- This imprecision together with the 
fact that (i?;)o is an approximation (see Eq. (05])) ex- 
plain the large absolute errors in Table HVl for parameters 
of high order. We shall detail the effects of the approxi- 
mation of (Ri)q' later in the charged case. 

As we have illustrated for the n-p system the numerical 
imprecision comes from big numbers given by R-matrix 
derivatives (see the previous value of r 4 ). This is not a 
general rule. Sometimes these kind of numbers comes 
from {"Hiijo ma king large values for t\ or t-i- There are 
three features to get these big numbers for the neutral 
case (see Eqs. (f2"Tj) and (|2"2"10 : 1) large channel radius, 2) 
big reduced mass and 3) high order of the effective-range 
parameter. 

When numerical inaccuracies appear for the effective- 
range parameters of high order, one can give a gross ap- 
proximation of these parameters. Let us illustrate it us- 
ing Table HVl For the last four effective-range parameters 
the imprecision decreases for small channel radii. It is 
totally expected following the argument number 1 pre- 
sented previously On the other hand, the precision for 
the first six parameters decreases for small channel radii 
(e.g., TV e {30,40} and a < 17 fm), which is due to 

I 



the approximation for the potential at r > a and not to 
numerical limitations. These two behaviors allow us to 
predict the order of magnitude for the first nine or ten 
effective-range parameters. It is achieved by choosing a 
small channel radius to avoid numerical inaccuracies for 
parameters of high order, but not too small to keep the 
correct order of magnitude for those of low orders. If an 
effective-range parameter is zero its order of magnitude 
is indeterminate, and therefore, we suggest to report the 
absolute error (e.g., for a = 5 fm and N = 40, the ab- 
solute error falls in the range [0.007, 0.7] fm n ~ for the 
first nine effective-range parameters of the n-p system) . 



B. Charged cases 



Comparing Fig. Q] and Fig. |4j one expects that the 
channel radius for the 3 He+a, 16 0+p and 12 C+a sys- 
tems will be larger than the channel radius of the W-S 
system. This intuitive prediction has been corroborated 
in Ref. for the 3 He+a and 16 0+p systems, where 
varying the channel radius from 10 fm up to 14 fm, the 
scattering lengths, the effective ranges and the shape pa- 
rameters are stable. Therefore, starting with the values 
of N and a which give us the best prediction according 
to Ref. [3|, and using the method described in Sec. Ill BL 
we display in Table |V] the numerical results for the first 
ten effective-range parameters of the systems 3 He+a and 
16 0+p. 



Table V. Results for the first ten effective-range parameters for 3 He+a and ie O+p. Units are shown in Table [ 









3 He+a (s-wave) 










O+p (p-wave) 






a (fm) 




12 




14 






14 




16 




TV 


30 




40 30 




40 


30 




40 30 




40 



111 

ri 
Pi 

Ql,3 

QlA 

Qi,& 
Qi,6 

Ql,7 

Qifi 

Ql,9 



36.8849 

0.97262 

-0.09010 

0.04972 

0.09278 

-0.07573 

-3.12 x 10" 

6.92 x 10 1 

6.21 x 10 

-1.65 x 



(i 
10 10 



36.8863 

0.97263 

-0.09009 

0.04970 

0.09281 

-0.07571 

-3.12 x 10" 1 

6.92 x 10 1 

6.21 x 10 6 

-1.65 x 10 10 



36.8847 

0.97262 

-0.09009 

0.04972 

0.09281 

-0.07581 

-7.45 x 10~ 2 

-4.93 x 10 2 

3.39 x 10 6 

-3.41 x 10 8 



36.8862 

0.97262 

-0.09008 

0.04971 

0.09282 

-0.075771 

-7.46 x 10" 

-4.93 x 10 2 

3.39 x 10 6 

-3.41 x 10 8 



401.99 
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-1.20 x 10 6 
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For both systems, the estimation of a/, n and Pi are 
in good agreement with the results of Ref. J9( . The next 
three parameter are stable and their values agree with 
the description of the phase-shifts at low energies, which 
can be checked by using Eqs. |7J and (J9)). These facts 
indicate that Qi : i is correctly estimated for i < 6. The 
last three parameters are evidently greater than those of 
lower order, which indicates that numerical limitations 
are present by using double precision calculations. 

The 3 He+Q; and 16 0+p systems do not have weakly 
bound states (our potential models provide bound states 



I 

at -34.2 MeV and -5.7 MeV for 3 He+a and -3.6 MeV 
for 16 0+p) and the phase-shits at low energies do not 
show indication about possible resonances. It is coherent 
with the stability of a; in Tabic [V] and corroborates that 
the imprecisions come from the numerical limitations. 

In order to detail the numerical limitations, we should 
track where the large numbers come from in the numeri- 
cal calculation. To do that we need to evaluate the contri- 
bution of the R-matrix derivatives and the renormalizcd 
Coulomb functions, in a similar way as we have shown 
in Eq. f4"5|) . Note that this equation is not valid for the 
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charged case, but removing the term n! from the denom- 
inator on the right side, one gets (-D; -1 )^ (see Eq. (|15jl). 
which has the same structure for the neutral and charged 
case. Thus, we can compute the terms ti, t%, t% and r 4 
and see if numerical imprecisions appear for determining 

For the 3 He+a system one finds (a = 14 fm, TV = 40) 



ti « -0.2 fm/MeV s 



5 x 10 11 fm/MeV a , 



T 2 

r 3 r -5 x 10 6 fm/McV 9 and r 4 ps 6 x 10 6 fm/MeV 9 . 
For the 16 0+p system one finds (a = 16 fm, TV = 40) 
ri ra 10 3 fm/MeV 9 , t 2 «4x 10 8 fm/MeV 9 , 
t 3 ps -7 x 10 3 fm/MeV 9 and t 4 ps 9 x 10 3 fm/MeV 9 . 

In Fig. [6] we present numerical results of (-Rz)o f° r 
j € {3,6,9}. This figure shows that the contribution 

of (Ri)q to form big numbers increases with j, as one 
could expect. However, despite rather large values of 
the R-matrix derivatives at zero energy ({Ri)q < 10 7 
MeV _J- if j < 9 and a < 16 fm for both systems), the 
values of T3 and r 4 indicate that there is not a relevant 
influence of these derivatives for numerical inaccuracies. 



Of course, this argument is valid if (Ri)q 
precisely (we shall discuss it later). 



U> 
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Figure 6. Numerical results of (Ri) for 3 He+a (square 
points) and ie O+p (circle points). Black (white) points cor- 
respond to TV = 40 (TV = 30). 

In contrast, Table IVII shows that the values of the 
renormalized Coulomb functions and their derivatives 
can be large enough to create numerical inaccuracies as 
seen by the large values of T2 with respect to those of t% 
and T4. Note that the results for n, T2, T3 and r 4 can 
be obtained partially by using Fig. |6] and Table IVII (the 
numbers in Table fVTl are rounded). 

Following Table [Vj numerical inaccuracies appear cal- 
culating Qo,6 f° r the 3 Hc+a system and calculating Qu 
for the 16 0+p system, i.e., to higher order for the lat- 
ter. It is not surprising from the numerical point of view 



Table VI. Values for the renormalized Coulomb functions and 
some of their derivatives at zero energy for 3 He+a, 16 0+p 
and 12 C+q. Ui\ n) is given in fm 1/2 /MeV". 





3 He+a 


le O+p 


12 C+a 




(s-wave) 


(p-wave) 


(d-wave) 


a (fm) 


14 


16 


17 


a N (fm) 


4.22 


3.83 


0.81 


iHu)^ 


186 


170 


1.3 x 10 5 


(*rf 


0.02 


0.02 


1.4 x 10~ 5 


{nut ] 


-299 


108 


-1.8 x 10 5 


{-Hi^ 


1.68 


0.1 


1.6 x 10~ 4 


(Hu)^ 


83.6 


-1016 


-5.8 x 10 5 


CHn)P 


9698 


81.6 


1.1 x 10 -2 


(nu)^ 


-7.4 


6.6 x 10 4 


5.9 x 10 6 


m 2l y* 


2.8 x 10 9 


2.1 x 10 6 


4.9 



because the numbers involved in the calculation of the 
cffcctivc-rangc parameters are smaller for the 16 0+p sys- 
tem. We think that this could be part of the explanation 
about why the parameter Qi,6 seems more stable than 
the parameter Qo,6 in Table IVl 

It is clear that the numerical inaccuracies are stronger 
for the 3 He+a system than for the 16 0+p system, and 
if the R-matrix derivatives are well computed, it is also 
clear that the most significant contribution to make big 
numbers comes from the renormalized Coulomb functions 
for both systems. Let us suppose for a moment that 
the numerical calculations were performed by using sin- 
gle precision. In this case, single precision cannot handle 
correctly additions and subtractions of T\ , T2, T3 and r 4 
for the 3 He+a and 16 0+p systems, and then we could 
stop the discussion at this point. However, all results 
shown here are obtained by using double precision, which 
can handle correctly those additions and subtractions. 
Therefore, the only possibility to get numerical impreci- 
sion comes from the approximation of (Ri)q , which can 
also be understood as a numerical limitation. 

As Eq. (|4"8]) stresses, a better estimation of the R- 
matrix derivatives at zero energy is achieved by increas- 
ing the Lagrange mesh points, TV. It means that there 
exists a value of TV, TV, for which one gets a "stable" 
approximation of (Ri) . Here, "stable" means that for 

TV > TV the value of (Ri) docs not change up to a 
specific digit in its fractional part. 

Let us illustrate the role that the previous effect plays 
for the 16 0+p system. Figure [6] shows that there is not 
a significant difference in the estimation of (Ri)q by us- 
ing TV = 30 or TV = 40. For these Lagrange meshes and 



16 fm, we obtain (i?i) 
\(9) 



(0) 



= 656.590 ■ • • MeV J and 
(Ri)q J> — 656.565 ■■ • MeV -9 respectively, which means 

that we have a precise determination of (Ri) up to 
the first digit in the fractional part. Our internal calcu- 
lations have shown that a precise determination of the 
second digit in the fractional part of (Ri) is possible 
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for N > 120. This effect makes a clear numerical limita- 
tion to determine the effective-range parameters by using 
the technique in Sec. IIII[ where in order to compute cor- 
rectly a large set of digits for (Ri)q, one should imple- 
ment high-precision in the numerical calculation because 
of the large mesh required. 

It is important to stress that the previous numerical 
limitation is remarkably strong when the value of (Ri) 
are large, which is often for determining effective-range 
parameters of high order (see n-p system as an example) . 

On the other hand, in the charged case there is an 
extra parameter, which is the nuclear Bohr radius a^ 
or the nuclear Rydberg energy Epf. Remember that En 
can be written in terms of aw (see Eq. ([3])). Hence, in 
order to explore the impact of a^, we have introduced 
the 12 C + a system, which has a nuclear Bohr radius 
five times smaller (approximately) than the nuclear Bohr 
radius of 3 Hc+a or 16 0+p. 

Table IVIII shows that the first ten effective-range pa- 
rameters associated to 12 C + a arc stable, which agrees 
with the results presented in Table IVII and the descrip- 
tion of the phase shifts following Eqs. (J7J and © at low 
energies. This stability and description are also expected 
from Eqs. (f3"5"j) and P0|) . where one sees that derivatives 
of the renormalized Coulomb functions are inversely pro- 
portional to some power of the nuclear Rydberg energy 
(see for instance Eq. (|36[) . which is necessary to com- 
pute Eqs. (|55|) and PUIO or directly proportional to some 
power of the nuclear Bohr radius. It leads us to infer that 
the numerical limitations introduced by the renormalized 
Coulomb functions will be more important when ajy is 
large or when En is low. 



Table VII. Results for the first ten parameters in the ERE for 
12 C+<x Units are shown in Table [I] 



(9) 



a (fm) 


15 




17 




N 


30 


40 


30 


40 


a 2 


58926 


58903 


58899 


58903 


r> 


0.15797 


0.15797 


0.15797 


0.15797 


Pi 


-65.964 


-65.964 


-65.964 


-65.964 


Q2,3 


0.14105 


0.14105 


0.14105 


0.14105 


Q2,i 


-0.00986 


-0.00986 


-0.00987 


-0.00986 


Q2,5 


-0.02996 


-0.02997 


-0.02997 


-0.02997 


Q2,S 


-0.02971 


-0.02973 


-0.02975 


-0.02973 


Q2,7 


-0.00154 


-0.00155 


-0.00157 


-0.00158 


Q2,S 


0.05553 


0.05557 


0.05550 


0.05547 


Q2,9 


0.13318 


0.13393 


0.13035 


0.12977 



Finally, we wish emphasize that the R-matrix and 
its derivatives for 12 C+a do not have large values to 
introduce strong numerical inaccuracies. For instance, 
for a = 17 fm and N = 40, (Ri)^ « 0.1 McV~ 3 , 
(i?;)o 6) « -5 McV~ 6 and (R)^ w 3 x 10 5 McV~ 9 . 
This means that for computing (£>/ _1 )o 
numbers such as n ~ 85 M f ^ , t 2 ~ 
r 3 « -4 x 10 6 ^ and r 4 



there will be 
6 x 10 MiV*' 



4 x lO 6 ^^, which can 



be checked partially, by using Table [VT1 and (Ri)q 



V. CONCLUSIONS AND PERSPECTIVES 

Calculations of the effective range parameters for neu- 
tral and charged nuclear systems are accurately carried 
out by solving the Schrddinger equation at zero energy 
for a given interaction potential, using the R-matrix and 
Lagrange-mesh methods. The expressions developed here 
show an easy path to go beyond the second order in the 
ERE without any mathematical or physical limitation. 
In particular, the present results are much more rigorous 
and systematic than methods based on a fit of low-energy 
scattering phase shifts, which typically lead to large error 
bars for high-order parameters 0] . 

Nevertheless, problems in the numerical calculation of 
the present method can arise due to computational im- 
precisions, in particular for large values of the channel ra- 
dius a. In general, choosing the minimal a such that the 
short range potential is negligible at distances larger than 
a, leads to a satisfactory gross estimate of the effective- 
range parameters. However, to improve the accuracy 
on the lowest-order parameters, larger radii are needed 
which lead to numerical instabilities for the high-order 
parameters. Using double numerical precision, the ex- 
pressions and methods explained here have worked very 
well up to the fifth or sixth order in energy on all exam- 
ples considered above. In some cases, a high accuracy 
could even be reached at least up to the ninth order. 

For typical meshes (30 to 40 mesh points) the 
Lagrange-mesh technique gives a fast convergence for 
the first digits of the R-matrix derivatives at zero en- 
ergy. However, the present method requires a large set 
of digits for these derivatives in order to compute the 
effective-range parameters precisely. This implies a much 
larger mesh which in general demands an implementation 
with high numerical precision. If such implementation is 
made, the numerical limitation given by large channel 
radius should be also decreased. 

When a weakly bound state is present a large channel 
radius can be demanded in order to keep the physical 
meaning of this state and the effective-range parameters. 
Similarly, if there is a virtual or resonant state at ener- 
gies very close to the threshold, choosing a large channel 
radius is mandatory in order to estimate the effective 
range parameters with a good precision. This case can 
thus be particularly delicate to handle since large radii 
can lead to numerical instabilities. However, when ap- 
plied to charged systems, our results have shown that 
the Coulomb barrier partly solves this problem. We infer 
that in general high-accuracy predictions of the effective- 
range parameters can be obtained without numerical dif- 
ficulty for nuclear systems with a large nuclear Rydberg 
energy or a small nuclear Bohr radius. The contribution 
of the centrifugal barrier could also play a secondary role 
in this improvement and would deserve a further study. 

The present work thus opens a window to get a bet- 
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ter description of elastic scattering for systems where the 
physical interest is at the low-energy regime, for instance 
in nuclear astrophysics. With this new tool, the problem 
of the connection between low-energy scattering phase 
shifts and weakly-bound-state asymptotic normalization 
constants (ANCs) could be revisited, at least for given 
interaction potentials. This is for instance important in 
systems of astrophysical interest like 12 C+a [7|. Once 
this connection is well understood for theoretical models, 
the possibility to directly extract ANCs from experimen- 
tal data could be considered. 

Another perspective is the use of Pade expansions 
rather than Maclaurin expansions for the effective-range 
function. These expansions are expected to be valid on a 
wider energy range and to allow the description of reso- 



nances, two useful features for a direct fit of experimental 
data [13|. These expected features could be checked on 
theoretical models first, hopefully leading to general pre- 
scriptions to prefer one type of expansion to another. 
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